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Abstract 

Cross sections for the production of two isolated muons up to high di-muon masses are 
measured in ep collisions at HERA with the HI detector in a data sample corresponding 
to an integrated luminosity of 71 pb" 1 at a centre of mass energy of y/s = 319 GeV. 
The results are in good agreement with Standard Model predictions, the dominant process 
being photon-photon interactions. Additional muons or electrons are searched for in events 
with two high transverse momentum muons using the full data sample corresponding to 
114 pb _1 , where data at y/s = 301 GeV and -y/s = 319 GeV are combined. Both the 
di-lepton sample and the tri-lepton sample agree well with the predictions. 
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1 Introduction 



Muon pair production in electron proton scattering proceeds mainly via two-photon interactions, 
77 — > where the incoming photons are radiated from the beam particles Q. It is impor- 

tant to check the quantitative agreement between experiment and theory in this process, since 
the understanding of this source of muons is vital in any search for anomalous muon produc- 
tion The clean experimental signature and the precise Standard Model prediction provide 
high sensitivity in such searches. In an analysis of multi-electron production @|, six outstand- 
ing events, three di-electron and three tri-electron events, were observed with di-electron masses 
above 100 GeV, a region in which the Standard Model prediction is low. A comparison with 
di-muon production in the same experiment is therefore particularly interesting. 

In this paper, a study of muon pair production in electron 1 proton scattering (ep — > efifiX) 
is presented using the HI detector at the ep collider HERA. The main part of this analysis is 
based on data with an integrated luminosity of 70.9 pb -1 collected with an electron energy of 
27.6 GeV and a proton energy E p = 920 GeV (y/s = 319 GeV). These data were recorded in 
the years 1999 and 2000 in positron proton scattering (60.8 pb -1 ) and electron proton scattering 
(10. 1 pb -1 ). Differential cross sections as functions of the invariant mass of the muon pair M m , 
the muon transverse momenta Pf 1 and the transverse momentum P* of the hadronic system X 
are measured for M m > 5 GeV. Results are also given for elastic and inelastic muon pair 
production separately. In addition, events with high di-lepton masses are studied in fifi and fifxe 
event samples with cuts adapted to the multi-electron analysis For this analysis the data at 
y/s = 301 GeV (E p = 820 GeV) and y/s = 319 GeV from the years 1994 - 2000 are combined, 
yielding a total luminosity of 1 13.7 pb -1 . 



2 Standard Model Processes 

The dominant process for the production of muon pairs in ep interactions is the two-photon 
reaction illustrated in the Feynman diagram shown in figure^. Due to the photon propagators, 
the momentum transfer to the scattered particles is generally small. While the calculation of 
the corresponding processes in e + e~ scattering is rather straightforward, the hadronic structure 
of the proton must be taken into account in ep scattering. The Feynman diagrams (b)-(e) in 
figure [l] show additional sources of muon pairs, which are less important than the two-photon 
reaction (a). They represent the four bremsstrahlung amplitudes in eq scattering with conver- 
sion of the radiated photon into a muon pair. The diagrams with the photons radiated from the 
electron lines can be viewed as Compton scattering of the electron and the photon, which is 
exchanged with the quark line (QED Compton process). In the phase space considered here, 
the Cabbibo-Parisi process, diagram (b), dominates (c) due to the pole of the Compton scatter- 
ing cross section for backward scattered photons in the ^e centre of mass system. A similar 
argument shows that the Drell-Yan process (d) dominates (e). Due to the negative interference 
of the diagrams (a) and (b) in the low mass region (M^ < 10 GeV), the expected cross section 
in the analysed phase space is about 5 % lower than that calculated when only the contribution 
of the two-photon process is considered. 

'in this paper "electron" refers to both electron and positron, if not otherwise stated. 
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Figure 1: Feynman diagrams for muon pair production in electromagnetic eq interactions, 
shown for the inelastic case where the scattering takes place from a single quark, such that 
the proton dissociates to a system X. In the elastic case, the proton scatters coherently. In 
the general electroweak case, any of the photon propagators can be replaced by a Z boson 
propagator. 



The program GRAPE [5| is used to compute the electroweak theory predictions. It calcu- 
lates the cross section according to the diagrams (a)-(e) of figure[T]utilising the GRACE [6] pro- 
gram. The program also includes contributions due to Z° boson exchange. These contributions 
become important in cases where the Z° radiated from external lines in diagrams (b)-(e) is close 
to mass shell. The w-pole contribution to the Drell-Yan-process (diagram d) is neglected. Its in- 
fluence was estimated in and was found to be negligible. Three different approaches are used 
to describe the proton structure in different phase space regions for the reaction ep — > e\i\iX. 
In the elastic region (X = p) the hadronic structure is parametrised by the electromagnetic 
form factors, which depend on Q 2 the negative four momentum transfer squared between 
the incoming and outgoing proton. In the quasi-elastic region (m p + < m x < 5 GeV) 
and the soft inelastic region (mx > 5 GeV and Q 2 p < 1 GeV 2 ), the calculation is based on 
parametrisations of the proton structure functions, which are given in for the nucleon reso- 
nance region (mx < 2 GeV) and in Q above the resonance region. In the deep inelastic region 
(m x > 5 GeV and Q 2 p > 1 GeV 2 ) the cross sections for eq — > e/jfiq are convoluted with the 
parton density functions of the proton lUUl. The GRAPE program is interfaced to PYTHIA [ 1 1 1 
and SOPHIA ITT21 for a complete simulation of the inelastic muon pair production processes. 
GRAPE allows for QED initial state radiation by adapting the cross section calculation in lfT3ll . 
Final state radiation is calculated with the parton shower method implemented in PYTHIA. 

Vector meson production with subsequent decay into muons is another source of muon pairs 
in ep scattering. Due to the mass cut (M w > 5 GeV) only the production of T mesons needs to 
be considered. The cross section is calculated using the Monte Carlo generator DIFFVM lfT4ll . 
Events with two muons also arise from the decay of tau-leptons produced in two-photon col- 
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lisions 77 — > tt and from semi-leptonic decays in open heavy quark production (QQ, i.e. 
cc and bb). These reactions are simulated with GRAPE and AROMA ifTSll . respectively. The 
leading order cross section for 6-production calculated with AROMA is normalised to the HI 
measurements presented in lfT*6l . 



3 Data Analysis 

3.1 Experimental Conditions 

The HI Detector ifTTl contains a central tracking detector (full acceptance over the range 
25° < 9 < 155°) and a forward 2 tracking detector (7° < 9 < 25°), which are surrounded by 
a liquid argon calorimeter (4° < 9 < 154°) and a lead-scintillating fibre calorimeter 3 (SpaCal 
calorimeter [18], 153° < 9 < 178°). The central tracking detector comprises proportional and 
drift chambers allowing a transverse momentum resolution of a(P t )/P^ = 0.005 GeV" 1 . These 
detector components are surrounded by a superconducting magnetic coil with a field strength 
of 1.15 T. The iron return yoke is equipped with streamer tubes forming the central muon de- 
tector (4° < 9 < 171°). In the forward direction, a proton remnant tagger (0.06° < 9 < 0.3°, 
z = 24 m) and a forward muon detector (3° < 9 < 17°) are used to separate elastic and inelas- 
tic processes. The trigger for this analysis is based on single muon signatures from the central 
muon detector, which are combined with signals from the central tracking detector. In events 
with large hadronic transverse momenta (P* > 12 GeV) trigger signals from the liquid argon 
calorimeter are used in addition. 

The procedure to extract cross sections relies on the HI detector simulation, which is based 
on the GEANT program JT9). After simulation, the generated events pass through the same 
reconstruction and analysis chain as the real data. Trigger and muon identification efficiencies 
are determined with high statistics data samples for the different subdetectors and are incor- 
porated in the simulation. Acceptances and trigger efficiencies for muon pair production are 
then determined from the Monte Carlo simulation and are used to correct the event yields to 
obtain cross sections. The overall trigger efficiency for di-muon events is about 70 %. At high 
hadronic transverse momenta {P* > 25 GeV) this efficiency is above 98 % 0. 

3.2 Event Selection 

The event selection and the muon identification are optimised in order to select an event sample 
consisting of two well identified muons, isolated from other objects in the event. The muon 
identification [20] is based on measurements from the central tracking detector, the central muon 
detector and the liquid argon calorimeter. Muon candidates are selected from tracks measured in 
the central tracking detector, which are linked to tracks measured in the muon detector. Muons 
which do not reach the muon detector, or enter inefficient regions of the muon detector, can be 

2 The forward direction and the positive z-axis are given by the proton beam direction. Polar angles 9 are 
defined with respect to the positive z-axis. The pseudorapidity is given by ?/ = — log(tan 0/2). 
3 This device was installed in 1995, replacing a lead-scintillator "sandwich" calorimeter 1171 . 
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identified by a central track linked to a signature of a minimal ionising particle in the liquid 
argon calorimeter. In about 10 % of the selected di-muon events, one of the muons is identified 
only in the calorimeter. The efficiency for identifying a single muon is typically 75 % in the 
kinematic range specified below. The momentum and angle measurements are obtained from 
the central tracking detector. 

The analysis requires two muons in the phase space given by: 

• polar angle region 20° < 9^ < 160°; 

• transverse momenta Pf 1 > 2 GeV and P/ 12 > 1.75 GeV; 

• invariant mass of the muon pair > 5 GeV. 

The polar angular range is matched to the acceptance of the central tracking detector, allowing 
for a precise momentum measurement. The requirement of a minimum transverse momentum 
ensures good muon identification. The analysis is focused on invariant masses above the J ftp- 
mass. Low muon pair masses are studied in H21U22I . 

Background from cosmic ray muons is suppressed by requiring that: 

• the ^-coordinate of the event vertex is within 40 cm of the nominal interaction point; 

• the opening angle between the two muons is smaller than 165°; 

• the timing of the event determined in the central tracking detector coincides with that of 
the ep bunch crossing; 

• the timing of the two muon track candidates is consistent with their emergence from a 
common vertex. 

The remaining cosmic background contribution is determined to be below 1 %. 

An isolation requirement suppresses events with muons from heavy quark decays and events 
with particles misidentified as muons: 

• the distance of the muons to the nearest track or jet 4 in the pseudorapidity-azimuth plane, 
Dj, rackJet = ^/ Arj 2 + A0 2 , is required to be greater than 1. Since at high transverse 
momenta the background contributions are small, the cut is relaxed to D^ rack Jet > 0.5 
for Pf > 10 GeV. 

The remaining contribution of misidentified muons is about 0.5 events, estimated using Monte 
Carlo simulations of photoproduction and neutral current deep inelastic scattering processes. 

The transverse hadronic momentum is measured ll2"4l using the liquid argon and SpaCal 
calorimeters, excluding energy deposits of identified muons or electrons. Electrons are identi- 
fied in the liquid argon or in the SpaCal calorimeter. 

4 Jets are considered with minimum transverse momentum of P/ et > 5 GeV, identified with the k t - 
algorithm 1231 . 
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3.3 Systematic Uncertainties 



The following uncertainties on the measured cross sections are taken into account. 

• The uncertainty on the integrated luminosity measurement is 1.5 %. 

• The uncertainty on the trigger efficiency, obtained from an independently triggered event 
sample, gives a contribution to the systematic uncertainty of 6 %. 

• The uncertainty on the identification efficiency of muons is determined by detailed com- 
parison of data and simulation efficiencies for a data sample consisting of events with 
exactly two tracks and at least one identified muon. This leads to a contribution to the 
systematic uncertainty of 6 %. 

• The uncertainty due to the reconstruction efficiency of the central tracking detector for 
the two muon tracks contributes with 4 %. 

• The uncertainties on the muon 9 and measurements are 3 mrad and 1 mrad, respectively, 
leading to an effect of up to 1 % on the cross section. 

• The systematic uncertainty due to biases in the transverse momentum measurement for 
high momentum tracks (P t > 20 GeV) is evaluated using the electrons in a neutral current 
deep inelastic scattering event sample. A scale uncertainty is derived from the ratio of the 
electron energy measured in the calorimeter to its track momentum measurement. The 
largest effect on the cross section is 7 % in the highest P t bin. 

These uncertainties added in quadrature lead to a total systematic error of 10 % on the integrated 
cross section. The uncertainty on the hadronic energy scale (4 % for the liquid argon calorimeter 
and 7 % for the SpaCal calorimeter) contributes an additional systematic error to the da/dP* 
determination. 

The uncertainty on the GRAPE calculation is below 1 % for the elastic process E51 . The 
accuracy of the calculation for the inelastic process is limited by the knowledge of the proton 
structure. The uncertainty on the structure function parametrisation (quasi-elastic) and the par- 
ton density function (deep inelastic) cause an uncertainty smaller than 5 %. The uncertainty on 
the predictions for other sources of muon pair production is estimated to be 30 % for QQ — > jiji 
and 40 % for T — > fi^i. For the error on the predicted event yields, these theoretical and the 
experimental uncertainties are added in quadrature with the statistical uncertainty on the Monte 
Carlo calculation. 

4 Results 

4.1 Inclusive Two Muon Cross Sections 

The cross section for the production of events with at least two muons is measured 5 and com- 
pared with the Standard Model prediction. In total, 1206 data events with two muons are se- 
lected, while 1197 ± 124 events are expected according to the Standard Model calculation. 

5 This analysis is based on data taken at a centre of mass energy of y/s = 319 GeV. More details can be found 
in ll20l . An analysis of data at ^fs = 301 GeV can be found in l26l . 
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No event with more than two muons is observed. In table [T] the contributions of the differ- 
ent Standard Model processes are given, where the errors contain the experimental and model 
uncertainties. The electroweak muon pair production process dominates all other processes, 
from which only 28.3 ± 6.7 events are expected. Of the 1206 data events, five events, all with 
M m < 11 GeV, have two equally charged muons. This observation is in good agreement with 
the expectation of 4.1 ± 1.5 events from the decay of heavy quarks. 

The cross section, evaluated in the phase space defined by > 5 GeV, Pf 1 > 2 GeV, 
Pf 2 > 1.75 GeV and 20° < 0^ < 160°, is presented in figureEK and tablets a function of 
the di-muon mass M w . The cross section falls steeply over more than four decades over the 
measured mass range, which extends up to 100 GeV. The shaded histograms show the expected 
contributions from the T and Z° resonances, where the latter is also included in the electroweak 
GRAPE prediction. At small masses minor contributions from open heavy flavour quark pro- 
duction, which are strongly suppressed due to the isolation requirement, and tau-decays are 
expected. The muon production cross section as a function of the transverse momenta of the 
two muons is presented in figure |2b an d table Both measured cross sections are in good 
agreement with the Standard Model expectations. The differential cross section as a function of 
the hadronic transverse momentum P* is also well described, as shown in figure|3]and tablelU 

The integrated cross section for electroweak muon pair production, a^w ^ | s obtained by 
subtracting the expected contributions from T, QQ and rr decays. The result is: 

= (46.4 ± 1.3 ±4.5) pb. 

The first error gives the statistical uncertainty and the second the systematic uncertainty. The 
measurement is in good agreement with the GRAPE prediction of (46.1 ± 1.4) pb. 

4.2 Elastic and Inelastic Muon Pair Production 

Elastic (ep — > e/i/xp) and inelastic (ep — > e/i/iX) muon pair production processes are distin- 
guished by tagging hadronic activity. An event is assigned as inelastic if activity is detected 
in the proton remnant tagger, the forward muon detector, or in the forward region of the liquid 
argon calorimeter (9 < 10°) [22|. Events containing tracks in the central or forward tracking 
detectors not associated to the muons or an identified electron are also considered as inelastic. 
A total of 631 data events are classified as elastic and 575 as inelastic. This is consistent with 
the Standard Model expectation, where 611 ± 87 elastic events and 586 ± 96 inelastic events 
are predicted. The Monte Carlo simulation shows that 92 % of generated inelastic events cause 
activity in the forward detectors and 93 % of generated elastic events remain untagged. The 
systematic uncertainty on the separation between elastic and inelastic pair production takes into 
account the tagging efficiencies of the forward detectors. This leads to an additional uncertainty 
of 10 % on the elastic and of 12 % on the inelastic cross sections. 

Figure Eland table |5] show the cross sections for elastic and inelastic muon pair production 
after subtraction of T, QQ and rr contributions. The two spectra are similar and are well 
described by the electroweak predictions. Elastic muon production contributes somewhat more 
in the low mass range and inelastic muon production has a slightly harder spectrum. This is 
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expected, as in elastic processes the electromagnetic form factors of the proton lead to a softer 
photon spectrum than that produced by radiation from point-like particles (inelastic process). 

In the analysed phase space an integrated cross section for elastic di-muon production of 

a% = (25.3 ± 1.0 ± 3.5) pb 
and for inelastic di-muon production of 

a™ 6 ' = (20.9 ± 0.9 ± 3.2) pb 

are measured. These measurements are in good agreement with the expected cross sections of 
(24.6 ± 0.3) pb and (21.5 ± 1.1) pb, respectively. 

4.3 Multi Lepton Events 

In addition to the determination of the inclusive di-muon cross section, events with two high P t 
muons and possible additional leptons, either muons or electrons, have been studied. In a small 
fraction of Standard Model electroweak di-muon production processes (figured), the electron 
is scattered through a large angle, such that it is visible in the detector and is not lost in the 
beam pipe, leading to an observed efifx final state. Events with three muons in the final state 
are suppressed within the Standard Model and a tri-muon signal would therefore be of great 
interest. In order to make use of the highest possible luminosity, data at y/s = 301 GeV are 
analysed in addition to the yfs = 319 GeV sample, resulting in a total luminosity of 113.7 pb -1 . 

To allow for a comparison with the multi-electron analysis [4>|, the following cuts are applied 
for this study: 

• two muons in the region 20° < 9 < 150°; 

• transverse momenta Pf 1 > 10 GeV and Pf 12 > 5 GeV. 

Additional muons must be detected in the central region of the detector, 20° < 6> M < 160°, 
with a minimum transverse momentum of 1.75 GeV. Additional electrons are searched for in 
the polar angle range 5° < 9 e < 175° and are required to have a minimum energy of 5 GeV. 
Such /i/xe events are triggered with an efficiency of typically 90 %. 

In the examined phase space, 56 di-muon events are found in the data, while 54.7 ± 5.7 
events are expected. Among these 56 events, 40 events contain exactly two muons (/i/i events), 
compared with 39.9 ± 4.2 expected. In the other 16 events (/i/ze events), one additional elec- 
tron is observed in the liquid argon or the SpaCal calorimeter, compared with an expectation of 
14.9 ± 1.6 events. As expected from the dominant two-photon process, the electron is preferen- 
tially found in the backward region of the detector. No event with three or more muons or with 
two muons and more than one electron is observed. 



10 



In figure |5k the di-muon mass distributions of events classified as /i/i events or as fifj,e 
events are compared with the theoretical expectations. Both invariant mass distributions are in 
agreement with the Standard Model calculations. The distribution in M 12 , the invariant mass 
of the two leptons with the largest P t , is shown for the /i/ie sample in figure |3J). This mass 
combination is selected in order to ease comparison with the multi-electron analysis H, where 
the scattered electron cannot be identified uniquely. For approximately half of these events, 
the two leptons with the highest P t are the electron and a muon. For these events, the mass 
distribution M^ 2 e is also shown in figure |3j5. Both mass distributions are compatible with the 
Standard Model predictions. 

For masses M 12 > 100 GeV (> 80 GeV) one n(i event is found, while 0.08 ± 0.01 
(0.29 ± 0.03) are expected. This inelastic event with two well identified muons has a mass 
of M M(U = 102 ±11 GeV and was recorded at E p = 820 GeV. No event classified as /i/ze with 
M 12 > 100 GeV is observed. The prediction is 0.05 ± 0.01. These results at high di-lepton 
masses are in agreement with the Standard Model predictions. In view of the present limited 
statistics, they cannot be used to draw firm conclusions concerning the high mass excess ob- 
served in the multi-electron analysis 0. 

5 Conclusion 

Isolated muon pair production is analysed for di-muon invariant masses above 5 GeV. The 
inclusive, elastic and inelastic cross sections are measured. In addition, a jijie event sample 
is studied. In all cases, the predictions of the Standard Model are in good agreement with the 
observations up to the largest di-lepton masses observed. 
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Inclusive Two Muon Production 



Inclusive Two Muon Production 




Figure 2: (a) Cross section for the production of two muons in ep interactions as a function of 
the di-muon mass M w . (b) Muon production cross section as a function of the muon transverse 
momenta P[ L (two entries per event). The data are compared with Standard Model predictions. 
See text for the accepted phase space. The relative difference between the data and the sum of 
all Standard Model contributions is also shown (lower figures). The inner error bars represent 
the statistical errors, the outer error bars the statistical and systematic errors added in quadrature. 
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Inclusive Two Muon Production 
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Figure 3: Cross section for two muon production as a function of the hadronic transverse mo- 
mentum . For further details see figure |2j 
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Figure 4: Cross section for electroweak (see text) muon pair production as a function of the 
invariant mass M M/J for elastically produced muon pairs (a) and for inelastically produced muon 
pairs (b) compared with the electroweak (EW) prediction using the GRAPE generator. For the 
accepted phase space see text. The inner error bars represent the statistical errors, the outer error 
bars the statistical and systematic errors added in quadrature. 
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Multi Lepton Analysis 



Multi Lepton Analysis 
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Figure 5: a) Distributions of the invariant mass M MAt for jiji (points) and jijie (triangles) events, 
compared with the Standard Model predictions, b) /i/ie event distributions of the mass Af 12 of 
the two highest P t leptons (points), and of the mass M^ 2 e for events where the leptons with the 
highest P t are a muon and the electron (triangles). For clarity, the closed data points are shifted 
slightly to the left. The error bars represent the statistical errors. 



Data 


SM 


EW fi+fx- 


T — > /i/x 


TT > /i/i 


QQ — > /i/i 


1206 


1197 ±124 


1169 ± 122 


12.3 ±5.1 


4.5 ±0.6 


11.5 ±3.8 



Table 1 : The number of selected di-muon events compared with the Standard Model prediction 
(SM). The dominant electroweak contribution (EW) is determined using the GRAPE generator. 
The expectations for other contributions are also given. 
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range 




[GeV] 


[pb/GeV] 


5.0 


- 5.3 


13.3 ± 1.3 ± 1.3 


5.3 


- 5.7 


12.8 ± 1.2 ± 1.2 


5.7 


- 6.0 


12.5 ± 1.2 ± 1.2 


6.0 


- 6.5 


9.09 ±0.89 ±0.86 


6.5 


- 7.0 


8.85 ±0.84 ±0.84 


7.0 


- 7.6 


7.13 ±0.68 ±0.68 


7.6 


- 8.2 


6.48 ±0.63 ±0.62 


8.2 


- 8.9 


4.11 ±0.45 ±0.39 


8.9 


- 9.8 


4.01 ±0.42 ±0.38 


9.8 - 


10.7 


2.65 ±0.33 ±0.25 


10.7 


- 11.8 


1.54 ±0.24 ±0.15 


11.8 


- 12.9 


1.36 ±0.22 ±0.13 


12.9 


- 14.1 


0.96 ±0.19 ±0.09 


14.1 


- 15.4 


0.88 ±0.16 ±0.08 


15.4 


- 17.1 


0.57 ±0.11 ±0.05 


17.1 


- 19.1 


0.369 ±0.089 ±0.035 


19.1 


-21.6 


0.165 ±0.052 ±0.016 


21.6 


-26.0 


0.090 ±0.029 ±0.009 


26.0 


-31.0 


0.075 ±0.025 ±0.007 


31.0 


-40.0 


0.027 ±0.011 ±0.003 


40.0 


-53.0 


0.0093 ±0.0054 ±0.0009 


53.0 


-70.0 


0.0019 ±0.0019 ±0.0002 


70.0- 


110.0 


0.00075 ± 0.00075 ± 0.00009 



Table 2: Cross section for the production of two muons in ep interactions as a function of 
the di-muon mass M w . Di-muon events from T, r-pair and QQ-decays are included in the 
measurement. The first uncertainty is statistical and the second systematic. 
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P^ range 


da/dP£ 


[GeV] 


[pb/GeV] 


1.8-4.0 


30.7 ±0.8 ±2.9 


4.0 - 6.2 


8.18 ±0.37 ±0.78 


6.2-8.8 


2.26 ±0.18 ±0.22 


8.8-12.5 


0.580 ±0.073 ±0.055 


12.5 - 17.5 


0.169 ±0.034 ±0.016 


17.5 - 25.0 


0.036 ±0.012 ±0.003 


25.0 - 40.0 


0.0136 ±0.0052 ±0.0014 


40.0 - 60.0 


0.0015 ±0.0015 ±0.0002 



Table 3: Muon production cross section as a function of the muon transverse momenta P£ (two 
entries per event). The first uncertainty is statistical and the second systematic. 



P* range 
[GeV] 


da/dP$ 
[pb/GeV] 


0.0- 12.0 
12.0-25.0 
25.0 - 40.0 
40.0 - 80.0 


3.94 ±0.11 ±0.38 
0.0174 ±0.0071 ±0.0029 
0.0027 ±0.0027 ±0.0005 
0.00097 ± 0.00097 ± 0.00023 



Table 4: Cross section for two muon production as a function of the hadronic transverse mo- 
mentum P 4 X . The first uncertainty is statistical and the second systematic. 
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dcr eZ IdM 


J inel lAj, .]- 


[GeV] 


[pb/GeV] 


[pb/GeV] 


5.0-6.8 


6.85 ±0.40 ±0.95 


4.19 ±0.29 ±0.64 


6.8-8.7 


3.16 ±0.25 ±0.44 


2.91 ±0.24 ±0.45 


8.7-10.6 


1.48 ±0.17 ±0.21 


1.74 ±0.19 ±0.27 


10.6- 12.8 


0.78 ±0.13 ±0.11 


0.60 ±0.10 ±0.09 


12.8-15.9 


0.375 ±0.065 ±0.052 


0.525 ±0.084 ±0.080 


15.9-19.3 


0.210 ±0.051 ±0.029 


0.207 ±0.050 ±0.032 


19.3 - 23.9 


0.073 ±0.026 ±0.010 


0.060 ±0.023 ±0.009 


23.9 - 30.0 


0.041 ±0.017 ±0.006 


0.039 ±0.016 ±0.006 


30.0 - 40.0 


0.0108 ±0.0063 ±0.0015 


0.0165 ±0.0082 ±0.0025 


40.0 - 55.0 




0.0092 ±0.0053 ±0.0014 


55.0 - 90.0 


0.00082 ±0.00082 ±0.00013 


0.00091 ±0.00091 ±0.00017 



Table 5: Cross section for electroweak muon pair production as a function of the invariant mass 
for elastically produced muon pairs (second column) and inelastically produced muon 
pairs (third column). Muons from T, r-pair and QQ-decays are considered as background and 
the expected event yields from these processes are subtracted from the measured event numbers. 
The first uncertainty is statistical and the second systematic. 
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